Multifractality of Hamiltonians with power-law transfer terms 



E. Cuevafl 

Departamento de Fisica, Universidad de Murcia, E-30071 Murcia, 
(Dated: February 2, 2008) 



Spain. 



o 
o 

(N 
> 

o 



I 



c 

o 

(N 
> 

(N 
O 

\o 
o 
m 
o 

-I— > 



T3 

c 

o 
o 



X 



Finite-size effects in the generalized fractal dimensions d 5 are investigated numerically. We con- 
centrate on a one-dimensional disordered model with long-range random hopping amplitudes in both 
the strong- and the weak-coupling regime. At the macroscopic limit, a linear dependence of d q on 
q is found in both regimes for values of q < 4g~ , where g is the coupling constant of the model. 
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I. INTRODUCTION 

Quantum phase transitions in disordered electronic 
systems remain one of the central problems in condensed- 
matter physics. Considerable attention has recently fo- 
cused on the critical eigenfunctions, which strongly fluc- 
tuate near the critical point and thus have multifractal 
scaling propertiesiS^i^ Wave-function statistics can be 
characterized through the set of generalized fractal di- 
mensions d q which are associated with the scaling of the 
gth moment of the wave-function intensity. A complete 
knowledge of d q is equivalent to a complete physical char- 
acterization of the fractali 

Metal-insulator transitions (MIT's) depend on the di- 
mensionality and symmetries of the system and can oc- 
cur in both the strong disorder and the weak disorder 
regime (strong-coupling or weak-coupling regime, respec- 
tively, of the corresponding field-theoretical description) . 
Each regime is characterized by its respective coupling 
strength g, which depends on the ratio between diagonal 
disorder and the off-diagonal transition matrix elements 
of the Hamiltonian. 6 

For a multifractal structure of wave functions, the q 
dependence of the fractal dimensions d q is given by the 
linear relation 



d q = d 



"gQ > 1 ^ b , 



(1) 



where c g <C 1 and b are parameters related to the cou- 
pling constant of the corresponding model (see below). 
This equation has been obtained for different transitions 
in the weak-coupling limit using different approaches, 
which we summarize in the following paragraphs. Equa- 
tion Q describes weak multifractality in the sense that 
the fractal dimensions d q slightly deviate from the sys- 
tem dimensionality d characteristic of the homogeneously 
spread metallic wave functions. 

In the context of the theory of Anderson localization, 
Eq. was obtained^ in the leading order in 1/2-kuD <C 
1, for two-dimensional (2D) conductors on the basis of a 
reduced version of the supersymmetric a model. In this 
case c g = \/2-k 2 v(3D, where v is the density of states, D 
the diffusion constant, and the symmetry parameter. 

Using a perturbative e expansion in one-loop order for 
a d = 2+e dimensional system, Eq. (JTJ was found^ in the 
leading order in e <C 1, with c g = e. This derivation was 



made for the case of unbroken time-reversal symmetry 
(orthogonal universality class) in the weak-coupling limit. 

Reference 8 derived Eq. with c g = g/8ir, in the 
limit g <C 1 for the power-law random banded matrix 
model (orthogonal universality class). This solution was 
found by mapping the problem onto a nonlinear a model 
with nonlocal interaction. 

Equation was also obtained, using the replica 
trick^ for the critical wave function of a 2D Dirac fermion 
in a random vector potential (chiral universality class). 
The c g = Ayi/V found depends on the variance of 
the Gaussian distributed vector potential. This result 
was confirmed recentljsifiiii by mapping the problem to 
a Gaussian field theory in an ultrametric space. 

Field-theoretical approaches in the context of the inte- 
ger quantum Hall transition (unitary universality class) 
have been addressed recently. Based on the application 
of conformal field theory to the critical point of a 2D 
Euclidean field theory, Bhaseen et a/. 1 — derived Eq. (JTJ 
with c g = 2/k, where k is the inverse coupling constant. 

We wish to point out that the existing theories are 
inherently weak-coupling results. However, MIT's gener- 
ically take place at strong disorder, where the energy 
scales associated with both disorder and interactions are 
comparable to the Fermi energy, unlike to what happens 
in perturbative 2 + e expansion approaches or in effective- 
field theories. Thus, quantitative results for d q are lack- 
ing at strong coupling, and finding d q for this regime is 
essential in order to fully understand the MIT's. 

The aim of this work was twofold. First, a numer- 
ical verification of the weak-coupling phase result, Eq. 
Q). Second, and most importantly, we were interested 
in finding the q dependence of d q in the strong-coupling 
regime, which has been left unexplored. Since this case 
is not accessible by the previous treatments, the problem 
can only be addressed by using numerical calculations. 

The paper begins by first giving the model and the 
methods used for the calculations in Sec. II. The results 
for the fractal dimensions in the weak- and the strong- 
coupling regimes are presented in Sec. Ill A and III B, 
respectively. Finally, Sec. IV summarizes our findings. 
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II. MODEL AND METHODS 

The disorder-induced MIT is usually investigated for 
Hamiltonians with short-range, off-diagonal matrix el- 
ements (e.g., the canonical Anderson model). Other 
Hamiltonians exhibiting a MIT in arbitrary dimension 
d are those that include long-range hopping terms^^ 
The effect of long-range hopping on localization was orig- 
inally considered by Anderson 1 ^ for randomly distributed 
impurities in d dimensions with the V(i — r') ~ \r — r'\~ a 
hopping interaction. It is knowniStik that all states are 
extended for a < d, whereas for a > d the states are 
localized. Thus, the MIT can be studied by varying the 
exponent a at fixed disorder strength. At the transition 
line a = d, a real-space renormalization group can be 
constructed for the distribution of couplings i 16 i 17 These 
models are the most convenient for studying critical prop- 
erties numerically since the exact quantum critical point 
is known (a = d) and, in addition, they allow the ID 
case to be treated, thus using larger system sizes and 
reducing the numerical effort. Many real systems of in- 
terest can be described by these Hamiltonians. Among 
such systems are optical phonons in disordered dielec- 
tric materials coupled by electric dipole forces^ excita- 
tions in two-level systems in glasses interacting via elas- 
tic strain^ magnetic impurities in metals coupled by an 
r~ 3 Ruderman-Kittel-Kasuya-Yodida interaction^ and 
impurity quasiparticle states in 2D disordered d-wave 
superconductors 

We will consider here the intensively stud- 
ied power-law random banded matrix model 
(PRBM) 4 i 5 i 8 i 22 i 23 i 24 i 25 i 26 i 27 . The corresponding Hamil- 
tonian, which describes a disordered ID sample with 
random long-range hopping, is represented by real 
symmetric matrices, whose entries are randomly drawn 
from a normal distribution with zero mean, (TH-ij) = 0, 
and a variance which depends on the distance between 
the lattice sites, 



<l^| 2 ) = r 



1 



* + 3 
i = 3 



(2) 



(l*-j|/&) 2 

where the parameter b plays the role of the inverse cou- 
pling constant b = 4g _1 of the corresponding a model 
at the center of the spectral band. The model de- 
scribes a whole family of critical theories parametrized by 
< b < co, which determines the critical dimensionless 
conductance, in the same way as the dimensionality la- 
bels the different Anderson transitions. In the two limit- 
ing cases b^> 1 and b <C 1, which correspond to the weak 
and the strong disorder limits, respectively, some critical 
properties have been derived analyticallyi^&SSiS&S 7 . More 
specifically, as already mentioned in the Introduction, in 
Ref. 8 the following equation in the limit 6^1 was 
derived, 



d q = l~— b , q «2nb 



(3) 



The system size ranges between L = 180 and L = 5400, 
and 0.01 < b < 12. We restrict ourselves to values of 



q < b, and consider a small energy window, containing 
about 5% of the states around the center of the band. 
The number of random realizations is such that the num- 
ber of critical states included for each L is roughly 4 x 10 5 , 
while, in order to reduce edge effects, periodic bound- 
ary conditions are included. Using methods based on 
level statistics, we checked that the normalized nearest 
level variances 2 ^ are indeed scale invariant at each critical 
point studied. 

For the computation of d q we used the standard box- 
counting procedure? 2 " first dividing the system of L sites 
into Ni = L/l boxes of linear size I and determining 
the box probability of the wave function in the i box 
by Pi(l) — I V'fcn | 2 ^ where the summation is restricted 
to sites within that box, and ipkn denotes the amplitude 
of an eigenstate with energy Ek at site n. The normal- 
ized gth moments of this probability constitute a mea- 
sure. From this, the mass exponents r q {L) 1 which en- 
code generalized dimensions d q (L) = r q (L)/(q — 1), can 
be obtained^ 



r q (L) 



lim 

.5^0 



(4) 



where S — l/L denotes the ratio of the box sizes and the 
system size. It should be made clear that the calculation 
of T q (L) is suitable only if the conditions 2 



a < I < L < £ 



(5) 



are satisfied, where £ is the localization or correlation 
length and a is the lattice spacing (or any microscopic 
length scale of the system). In practice, T q (L) is found 
by performing a linear regression of hiX^i-P? (0 with 
In 8 in a finite interval of 8 (usually a/L < 8 < 1/2). In 
order to properly satisfy the previous conditions JSJ and 
(£fj, in this work we will compute the derivative 



r q (L) 



d\nS 



(6) 



5=0.1 



First, T q (L) was calculated for different system sizes 
and then extrapolated to the macroscopic limit r q = 
liniL^oo T q (L). 



III. RESULTS 

A. Weak-coupling regime (b 3> 1) 

Using the exact eigenstates of Hamiltonian (J2J ob- 
tained from numerical diagonalizations, we evaluate, for 
each value of q, L, and b, the numerator on the right-hand 
side of Eq. for decreasing box sizes, and then calcu- 
late T q (L) from the slope of the graph of the numerator 
vs In 8 at 8 = 0.1, Eq. JJjJl. Fig. 1 provides an example of 
the In 8 dependence of ln^ i J 1 pj(l) in the weak-coupling 
regime (b — 10) for a system size L — 5400 and differ- 
ent values of q: 4 (circles), 5 (triangles), 6 (diamonds), 
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FIG. 1: The In 6 dependence of hi^i^iPiffl m the weak- 
coupling regime (6 = 10) for L = 5400 and different q val- 
ues: q — 4 (circles), 5 (triangles), 6 (diamonds), and 7 
(squares); the dotted vertical line corresponds to S = 0.1. 
The inset shows the same dependence of corresponding slopes, 

<nn££ipj(0/«nn«. 



and 7 (squares); the dotted vertical line corresponds to 
S = 0.1. The inset shows the same dependence of the 
corresponding slopes, rflnX^Ji Pi {l)/dlnS. Note that 
these derivatives are practically constant in the region 
shown. We have checked that the results were practi- 
cally the same when T q (L) was obtained from the linear 
fit to lnX^xPfCO vs In 5 in the interval 0.1 < S < 0.4. 

The size dependence of T q (L) at the critical point of 
the finite system for b = 10 is shown in Fig. 2 for q = 4. 
Clearly, significant finite-size effects are present, and for 
the larger values of L the convergence is evident. In order 
to predict the asymptotic value of r q from such a plot, 
a curve of the form r q (L) = r q + F a c (L) was fitted to 
the data points in this graph, where r 9 , a q , and c q are 
three fitting parameters and F aq , c „ (L) is a function which 
tends to zero as L — > oo. Various forms for F a iC (L) 
were chosen and tested on the data plots, including ex- 
ponential (~ e~ CqL ), inverse logarithmic (~ l/\nc q L), 
and power-law (~ L~ Cq ) decay with L. The testing in- 
volved performing a three-parameter fit on the T q (L) vs 
L data plots using the standard Levenberg-Marquardt 
method for nonlinear fits. The form of F aqtC (L) eventu- 
ally chosen was the power law. The reason for choosing 
the above form with preference to the alternative forms 
was simply that the inverse logarithmic and the expo- 
nential laws did not fit our data properly. In addition, 
for all values of q and b we found that the exponent c q 
hardly differs from unity. Thus, we can write 

T q (L) =r q +a q /L , (7) 

and keep only two free parameters. The solid line in Fig. 
2 is a fit to Eq. Q . Note that this equation gives a fairly 
good fit to the data. r q (L) is represented on a 1/L scale 
in the inset of Fig. 2 for different values of q: 4 (circles), 




L 

FIG. 2: The L dependence of the mass exponent r g (L) in the 
weak-coupling regime (b = 10) for q = 4 (circles); solid line 
is a fit to Eq. Q. The inset shows r q {L) on a 1/L scale for 
different q values: q — 4 (circles), 5 (triangles), 6 (diamonds), 
and 7 (squares); the straight lines are linear fits to Eq. Q. 



5 (triangles), 6 (diamonds), and 7 (squares). Here, one 
can better appreciate the linear behavior, Eq. Q , of the 
finite-size corrections to r q . 

The q dependence of the extrapolated results of d q = 
Tq/(q — 1) in the weak-coupling regime are depicted in 
the inset of Fig. 3 for several values of b: 6 (squares), 
8 (triangles), 10 (diamonds), and 12 (circles). The solid 
lines are fits to the form d q = 1 — c^q. The corresponding 
6-dependent slopes ct, are summarized and compared with 
the nonlinear er-model estimates l/27r6 in Fig. 5. The 
main panel shows the same data d q as a function of the 
rescaled variable q/2irb, in which it is clearly seen that the 
data points collapse fairly well onto a single straight line 
according to Eq. © (solid line). To our knowledge, this 
result provides the first numerical test of the nonlinear er- 
model prediction for the PRBM model. A similar linear- 
behavior has been found numerically for the distribution 
of the local electric fields at a dielectric resonance in two 
dimensions^ 



B. Strong-coupling regime (b -C 1) 

We have also calculated d q for the critical wave func- 
tions of Hamiltonian in the strong-coupling regime for 
which, as we already mentioned, there are no analytical 
predictions for values of q < b <C 1. The only analyt- 
ical estimate in this regime was obtained for q > 1/2, 
based on a resonance pairs approximation, in Ref. 22. 
r q ~ 4bT(q — ^)/y/nT(q — 1) found diverges at q = 1/2 
and for smaller values of q this approximation completely 
breaks down. 

Figure 4 shows the same extrapolated values of d q as 
in Fig. 3 in the strong-coupling regime (b <C 1). The 
values of b reported are 0.05 (squares), 0.1 (triangles), 
0.15 (diamonds), and 0.2 (circles), and the corresponding 
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q/2nb jj 



FIG. 3: d q as a function of the rescaled variable q/2nb in the 
weak-coupling regime (b 2> 1) for several lvalues: 6 (squares), 
8 (triangles), 10 (diamonds), and 12 (circles); the solid line 
represents the nonlinear a- model estimate, Eq. JJ. The inset 
shows the same data on the scale q; solid lines are fits to the 
form d q = 1 — Cbq. 




ql(A+Cb) 

FIG. 4: As in Fig. 3, for the strong-coupling regime (b <C 1): 
b = 0.05 (squares), 0.1 (triangles), 0.15 (diamonds), and 0.2 
(circles). The solid line corresponds to Eq. (|HJ. The straight 
lines in the inset are fits to the form d q — 1 — ctq. 



FIG. 5: The inverse slope l/c& (solid circles) of the fractal 
dimension as a function of the parameter b of the PRBM 
model. The solid lines represent the analytical asymptotic 
(b 3> 1) result 2nb and the strong-coupling (b <C 1) result 
proposed in this work, A + Cb. The small-6 behavior of 1/q, 
is reported in the inset; the straight line is a linear fit to the 
form A + Cb. 



paper. Thus, our results suggest the intriguing possibil- 
ity that the linear dependence of d q on q seems to be a 
generic law valid in all regimes. Although this is a conjec- 
ture at the present state, we expect that it will stimulate 
further analytical work. 

The inverse slope 1/q, of d q is depicted on a log- log plot 
in Fig. 5 (solid circles) as a function of parameter b of the 
PRBM model. The inset shows the small-6 dependence 
of 1/ch, from which the linear behavior on b is evident; 
the straight line is a linear fit to the proposed form with 
parameters A = 0.580 ± 0.009 and C = 4.17 ± 0.08. The 
numerical calculations agree perfectly with the analytical 
asymptotic in the weak-coupling regime 6 > 1, 1/c;, = 
2-7T&, and with the form proposed in this work for the 
strong-coupling regime b <C 1, 1/cfe = 0.58 + 4.176. 



IV. SUMMARY 



slopes Cb are summarized in Fig. 5 (solid circles). As in 
the weak disorder case, a linear dependence of d q on q 
is found for q < b. In order to collapse all data points 
onto a single straight line, as was done for the opposite 
limit b ^> 1, we propose the relation 1/q, = A+Cb for the 
slopes, where A ^ and C are two fitting parameters (see 
inset of Fig. 5). The collapse of the data is evident from 
the main panel of Fig. 4, when these data are represented 
as a function of the rescaled variable q/ (A + Cb). Hence, 
the q dependence of d q in the limit b <C 1 can be described 
by the linear relation 



dg = l~ 



A + Cb 



(8) 



This equation constitutes the main result of the present 



In this paper we have calculated the fractal dimen- 
sions d q of the wave functions for ID disordered systems 
with long-range transfer terms at criticality. The leading 
finite-size corrections to d q decay algebraically with ex- 
ponents equal to —1. At the infinite-size limit, we have 
confirmed that, according to theoretical predictions, d q 
is linearly dependent on q in the weak-coupling regime. 
Our calculations strongly suggest that this behavior is 
also valid at strong coupling. So, the linear dependence 
of d q on q seems to be a generic law valid in all disorder 
regimes. 

The question arises as to whether these results are ap- 
plicable to other quantum systems, particularly in the 
3D Anderson transition that occurs in the strong disor- 
der domain, and whose similarity with the PRBM model 
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at b = 0.3 has been demonstrated for several critical and BFM2003-03800. 
magnitudes^ 
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